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Abstract 

We approximate boson thermodynamic integrals as polynomials in two vari- 
ables chosen to give the correct limiting expansions and to smoothly interpolate 
into other regimes. With 10 free parameters, an accuracy of better than 0.009% is 
achieved for the pressure, internal energy density and the number density. We also 
revisit the fermion case, originally addressed by Eggleton, Faulkner and Flannery 
(1973), and substantially improve the accuracy of their fits. 
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1 INTRODUCTION 

The thermodynamic functions of fermions and bosons are expressed as in- 
tegrals which cannot be evaluated analytically except in limiting cases of 
degeneracy or relativity. Numerical integration is not entirely straightfor- 
ward, however, since the form of the integrand varies widely under different 
conditions. Thus the numerical method chosen must be carefully tested and 
optimized for the parameters at hand and different methods may be needed 
in different regimes. Furthermore, the computational time required can be 
an issue if the integrals need to be evaluated many times. For fermions, 
these problems can be overcome by using an ingenious numerical approach 
described by Eggleton, Faulkner, & Flannery (1973), hereafter referred to as 
EFF. They approximate the integrals by means of a polynomial in two vari- 
ables. The polynomial form is chosen so that it yields the correct behavior in 
the four limiting situations for which the series can be derived analytically. 
The coefficients are then optimized to accurately interpolate into regimes 
which are inaccessible to analytical analysis. We have found that a sim- 
ple modification of the approach of EFF, which involves a single additional 
parameter, can significantly improve the accuracy for a given number of co- 
efficients and, in addition, permits the approximation to merge exactly into 
the correct asymptotic behavior. This modification is superior to the one 
discussed by Pols et al. (1995). We discuss the fermion case briefiy in §3. 
The main purpose of the present paper is to analyze the thermodynamic 
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integrals for bosons; any condensate contribution would, of course, have to 
be handled separately. While we follow the general approach of EFF, bosons 
require a different approximation scheme. This is discussed in §4, in which 
we also present results for various orders of approximation. 

We begin by defining the integrals at hand in §2. The fcrmion and boson 
cases arc discussed in §3 and 4, respectively, and brief concluding remarks 
are given in §5. 



2 DEFINITIONS 



For a relativistic system of particles, the standard form for the pressure is 



±9T 



{27Thy 



In 



lie 



-(Efc-/.)/T 



(1) 



where g is a. degeneracy factor, the energy — ^JK^c^ + m^C^, is the 
chemical potential and the upper (lower) sign refers to fermions (bosons). 
The units employed are such that Boltzmann's constant ks — 1, so that 
the temperature is measured in units of energy. The angular integration is 
trivial for an infinite system and it is convenient to perform the standard 
integration by parts. Expressing the result in terms of the dimensionless 
variables of EFF 



t 



T 



fi — mc 



(2) 



and a dimensionless integration variable I — k/mc, we write the pressure in 
dimensionless form 



dl 



\/FTT 



(3) 
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where E — \JV- + 1 — 1 and d = ^mc^ (j^) carries the dimensions of 
P. For fermions the value of ip, or equivalently /i, is unrestricted. However 
for bosons, since the occupation probabihty has to be a positive quantity, 
ipt = — 1 < 0; thus ip is always negative. For a gas containing antibosons 
(b) in thermal equilibrium with bosons (b), so that /ij = — the requirement 
"06^ < yields the further restriction ipht = ipt > —2. 

The remaining thermodynamic functions are easily obtained. We write 
the entropy density S in dimensionless form as 
mc^S f dp\ 



s — 



d 



dt 



(4) 



The number density Af written in dimensionless form is 

dl f 



mc 



p 



d 



i = 7 



et-'^ ± 1 



(5) 



The energy density £ written in dimensionless form is 



e^ — = ts-p + piipt + 1) 
d 



dl l^Vf 



± 1 



-1 



(6) 



For e, the leading power series behavior of p obviously differs from that of 
ts, p and pipt. We therefore work with the internal energy density 



u 



e-p 



(7) 
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3 APPROXIMATION TO FERMION 
INTEGRALS 



3.1 Formalism 

The approach we use follows that of EFF. We determine an approximation 
for the pressure in terms of the degeneracy parameter and the temperature 
parameter t. Then, by eqs. (^H^), the remaining thermodynamic variables 
can be found by differentiation. The EFF scheme requires knowledge of the 
expansions of the pressure in the limiting cases. In the non-degenerate case 
(ND), ip —1, the denominator in equation (H) can be expanded in powers 
of exp^ip + 1/t) (Chandrasekhar 1958), leading to 



where, as before, the upper sign refers to fermions and the lower to bosons. 
Here K2 is a modified Bessel function which may be expanded (Abramowitz 
& Stegun 1965) in the extremely relativistic (ER, t ^ 1) and non-relativistic 
(NR, t <^ 1) limits. In the extremely degenerate case (ED) for fermions, 
in which ip ^ 1, the Sommerfeld expansion (Chandrasekhar 1958) may be 
employed and the extremely relativistic (ER, ipt ^ 1), and non-relativistic 
(NR, ipt <^ 1) limits obtained. Then the four limiting cases for the fermion 
pressure are 
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4^2, , ,5 f 45 , 15 , , Stt^ 1 

+ ^ + f + Pp + iS + l^ + -}.EDER (9) 

with the corresponding hmiting cases for the number density and internal 
energies easily obtained from equation by differentiation. 

The key to the EFF scheme is to find functions flip) and g{ip,t) such 
that equation can be expressed as 

oo oo 

P = fg'^' E E a^„/™^- ND NR ^7 « 1, / « 1, 

m=0 ?i=0 



oo oo 



P = f9' E E hmn^g-^ ND ER (7 » 1, / « 1, 

m=0 n=0 

00 00 

P = g"-'^ E E c„„/-'"(7" ED NR ^ < 1, / > 1, 

m=0 n=0 

00 00 

P = g^ E c^™„/-'"^7^" ED ER (7 » 1, / » 1, (10) 

m=0 n=0 

where Omn, &mn, and c/mn are coefficients. This is possible provided that 

00 00 
/(^) = e^ Ume"^^ ■ giij,t) = t E v^e"'^ ND < -1, 

m=0 m=0 

00 00 

/W = V-' E wm^-'™ ; = E x^V'-'" ED ^ > 1,(11) 

m=0 m=0 

where Um, Vm, Wm and Xm are additional coefficients. Truncating the sum- 
mations at order M and A^, equations (0) can be combined in the single 
expression 

which will be used to interpolate into all regions of and t. 
Equation (|TT]) indicates that / should obey 



/ cx e'^ ; cx /, ^ « -1, / « 1 ND, 

dip 
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/oc^2 ; ^«\//, /»1 ED. (13) 



We choose for the derivative 
df f 



(14) 



which differs from EFF by leaving free parameter to be fitted, rather 

than taking it to be unity as they did. This simple change gives a significant 
improvement in the fit, because it allows greater freedom in matching the 
limiting expression for the entropy in the extremely degenerate case. In 
the EFF scheme the relative errors in the entropy in the degenerate limit 
remain substantial even when the order of the approximation is very large. 
Integration yields 



^ = 2^ 1 + ^ + lnL^ . (15) 

« • /l + ^ + 1 



For the non-degenerate case, / ~ Aae"^'"^, which is small; while for extreme 
degeneracy, / ^ which is large. We define g = t-^/T+y, as in EFF, 

since we have found no advantage in introducing additional parameters here. 

The leading terms of equation (^ yield values for the four "corner" coeffi- 
cients and these are collected in Table 1. In addition, the leading contribution 
to the entropy density for the extremely degenerate case supplies constraints 
on the Pm-1,0 and Pm-i,n coefficients. These are also listed in Table 1. 

We choose the internal energy density and number density to be thermo- 
dynamically consistent with equation (1T^), so by differentiation we find 



J M N 

u = — — E Pmnrg" 

m=0 n=0 
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+ (3 _ l^A ll 



+(i + 1^ - A^)t^ + (i - 1^) ,^ , /J^ , ^ . (16) 



3.2 Results 

We determined the polynomial coefficients Pmn by a least squares fit to data 
for the pressure, which was obtained by accurate numerical integration using 
the same (/, g) grid as in EFF. This was done in order to facilitate compar- 
isons with EFF. The accuracies of the number density are of the same order 
as those for the pressure, so this procedure yields acceptable results for the 
internal energy and number density as well. Table 2 displays the maximum 
modulus (MM) and root mean square (RMS) deviations for various values of 
M and N. Specifically these are defined in terms of the fractional deviation 
of, for example, the pressure at the fitting points i, Apj, according to 

„Poly _ exact /VYAr).)2 

^Pi = ^exact' ; MM = Max| Ap,| ; RMS = y . (17) 



Table 2 also shows that treating Qj clS db free parameter significantly improves 
the fit as long as both M and N are greater than 1. For example, for the 
M — N — 3{2) case the improvement amounts to a factor of 8 (4). For 
fermions the number of coefficients to be fitted is {N + 1)(M + 1). 
The N = M = 2 case, whose coefficients are listed in Table 3, shows RMS 
deviations of 0.03%, which may be sufficient for many purposes. However, an 
order of magnitude improvement in the accuracy can be obtained by going 
to the N — M — 3 approximation, which is displayed in Table 4. 
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In many situations, it is preferable to have approximations that smoothly 
merge onto the exact limits. This prevents discontinuities from occurring 
when the exact results are used instead of the approximate scheme in extreme 
situations. We performed additional fits in which the 6 constraints from Table 
1 were incorporated. For the M = iV = 3 case (with 10+1 free parameters), 
the deviations, denoted by a dagger in Table 2, are about 4 times larger than 
in the unconstrained fit. Nevertheless, the overall accuracy is still acceptable 
for most numerical work, and if more accuracy is desired, one could increase 
M or N. Note that for the case M — N — 3, the improvement in the pressure 
error gained by fitting a amounts to a factor of 9. We list the coefficients for 
this case in Table 5. 

4 APPROXIMATION TO BOSON 
INTEGRALS 

4.1 Limiting Cases 

The approximation to boson integrals follows by analogy to that for fermions. 
However, the degenerate case is now denoted by ^ 1 since ip must always 
be negative. Furthermore, a boson-antiboson gas in equilibrium has the 
additional restriction that ipt > —2, which would preclude the existence of 
the non-degenerate, extremely relativistic case, although it is nevertheless 
useful to consider this formal limit below. 

We first determine the expansions in the various limiting cases. 
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4.1.1 Non-degenerate 

The expression for the pressure in the non-degenerate {ip <^ —1) case was 
given in equation and from this the internal energy density and density 
are easily obtained. 

For the non-relativistic case, in which t <C 1, the asymptotic series for 
the modified Bessel functions yields 




ND NR (18) 



As we have remarked, we consider the extremely relativistic (t ^ 1) case 
even though it is not physical when one is dealing with a boson-antiboson 
gas in thermal equilibrium. Making the small argument expansion of the 
modified Bessel functions, one obtains 
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4.1.2 Extremely Degenerate 

In the non-relativistic case {{ipl <C l,t <^ 1), we use y = E/t as the integra- 
tion variable in equation (j^) to find 

oo 

p = \te^ J dy {2yt)i {l + \yt + l^{ytf - ■ • •} [e^ - e^Y' . (20) 



The integrals can be written in terms of a function $ which is defined by 
Erdelyi et al. (1953) as 

oo 

$(e'^,s,l) =[V{s)]-^ j dtt'-\e' -e^]~^ ; s > 1, ^ < 



= |r(l - s){-^r' + E C(^ - r)^| , (21) 

I r=0 J 

where ({z) is the Riemann zeta function. For the series expansion to be valid, 
{ipl < 2-7? and s ^ 1,2,3. . .. This function is related to the polylogarithm 
function employed by Haber & Weldon (1982a, b); $(e^, s, 1) = e-^Li,(e'^). 
Using equations (pO| , |2TD we find 

P = (i) + C (1) ^ + K (i) + |0F(-^)^ + f C (i) t 

+fc(i)v^t+ic(i)t^ + ---} , 

" = (i) + ^ (i) ^ + K (I) + Iv^(-V^)^ + f C (i) t 

+fc(i)v^t+fic(i)t^ + ---} , 

P = (i) + C (^) + K (-1) - 2v^(-^)^ + f C (I) t 

+fC(i)V^t+TiC(i)t' + |v^t(-^)^ + ---} • EDNR (22) 

The opposite case of extremely relativistic, degenerate bosons {\ip\ -C 
l,t ^ 1) requires a series expansion of the pressure integral, followed by a 



12 



Mellin transformation and the evaluation of the residues at the poles, 
the results of Haber & Weldon (1982a, b), we find 

' -2 1 \ lla^ 7a 



Using 
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p = 2*« C(4)+C(3)a. + C(2) y-j^ + 36 
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where a = ^ + t'K Note that -1 < at < 1, if -2 < < 0. 
4.2 Numerical Approximation 

We now follow the strategy employed for fermions by EPF. The preceding 
formulae suggest that for the non-relativistic (relativistic) case, we need a 
power series in t {t~^). Thus, a suitable choice is simply g = t. In the non- 
degenerate limit, powers of are required. For the degenerate case, powers 
of are indicated by equation (^2]); the situation is less clear for equation 



(^31) , which has a complicated structure, but it would seem reasonable to use 



13 



a similar approach. We are therefore led to consider a variable h such that 

hoce-"^ ; ^ oc -h, 1^1 > 1, /i > 1 ND 
dip 

dh 1 



hoc^-ip ; ^oc-;^. |^/'|<1, /Kl ED (24) 

This does not uniquely specify the derivative and, after some experimenta- 
tion, we chose 

dh {^/a + hf 

# = h ' (2^) 

which has the desired limiting behavior. The introduction of the parameter 
a allows greatly improved fits, as we have also found in the fermion case (see 
§3). This property is related to the limiting behavior of the number density 
and entropy in the degenerate limit. Integration then yields 



^ = ^-ln ^ . (26) 

^+h \ ^ J ^ ^ 

For small values of h and /i ~ ^ —1ai\}^ while for large values h ~ ^Jae^~'^\ 
note that h > for all ip < 0. 

We then write the polynomial expansion for the pressure 

tf(l + t)f M N 

where the coefficients of the polynomial are denoted by Pmn- (Note that in 
the case h = 0, the interpretation h^ = 5om is to be made.) For the ED case, 
where h and are small, p should not contain a term proportional to h, 



which would yield a leading divergent l/y/—ip contribution to the number 
density. Therefore, we must impose the condition pin = (M + l)pon so that 
the contribution to equation (^) which is linear in h vanishes. Equation (p^) 
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then has the correct behavior in the four hmiting cases and should be suitable 
for smoothly interpolating for all values of ip and t. Since it is desirable 
that the internal energy density and number density be thermodynamically 
consistent with the pressure, we differentiate equation ( pTf ) to obtain 

tf(i + t)f _.„r, t 



7/ = - 



m=0 n=0 
i M N 



l + t 



P = + E ^^Pmnh-H-[-m + MM + 1 - m)]. (28) 

The entropy density can then be obtained from s = [p + u)/t — pif). We 
note that the condition pin = (M + l)pon leads to explicit cancellation of the 
leading order contributions to the entropy in the degenerate limit. 



Note that equation (|27D for bosons is actually the same as equation (|l^) 
for fermions as may be seen by setting f = 1/h and using the boson pa- 
rameterization g = t. It is simply more convenient to use these forms so 
that infinite values of h for bosons or / for fermions are only encountered for 
infinitely large In the boson case, it is necessary to be able to specify the 
case ip = 0, which is just h = 0. 

The leading terms for the limiting expressions for p and p, given in §4.1, 
yield values for selected coefficients Pmn and these are collected in Table 6. For 
reference, we note that C(3/2) = 2.61238, C(5/2) = 1.34149, C(3) = 1.20206, 
and C(4) = 7r^/90. The fitted values for these coefficients should approach 
the values in Table 6 with increasing dimensions M and N, as indeed we find. 
For the case M = 2 we have two different requirements on the coefficients, 
derived from the limiting cases for p and p. For p2o these expressions agree 
for a = 0.9486 or a = 0.9999, while for p2N they agree if a = 0.3854 or 
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a = 0.9005. Using a = 0.91, for example, the expressions agree to < 0.25% 
in both cases. On the other hand, for a = 1 the expressions for agree to 

0. 0003%, but the ones for p^o agree to only 2.7%. This large error explains 
why it is advantageous to employ a as an additional fitting parameter. 

4.3 Results 

The polynomial coefficients pmn in equations (p7| ) and ( pH]) were obtained 
by a least squares fit to "exact" results obtained by numerical Gaussian in- 
tegration (the accuracy was at least six decimal places). The results were 
generated for a (t, ip) grid consisting of t values of .05, .1, .5, 1, 5, 10 and 100 
and, for each t, seven equally spaced values of ip in the range —2<ipt< 0. 
The overall quality of the results obtained for various dimensionalities M 
and is indicated in Table 7 in terms of the maximum modulus and the 
root mean square deviations, as for fermions. If the coefficients pmn are ob- 
tained by fitting just the pressure data, then a comparable fit is obtained 
for the internal energy. However the fit for the number density can be sig- 
nificantly improved by separately optimising the pmn coefficients. Therefore 
the numbers in Table 7 were obtained by simultaneously fitting the pressure 
and number density data, thus sacrificing some accuracy in p to achieve im- 
provement in p (both the sacrifice and improvement are typically a factor of 
~ 4 in the RMS deviations). In addition to minimizing with respect to the 
(A^ + 1)M coefficients Pmn, the results were minimized with respect to the 
parameter a of equation (|26|) . Although a turned out to be fairly close to 

1, there is a substantial improvement gained by allowing it to vary. For the 
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cases listed a was in the range 0.7 — 1.1. 

The simplest approximation that is sensible, M — 2 N — 1, produces 
results accurate to better than 5% with only 4+1 parameters, while the M — 
N — 2 case with 6 + 1 parameters achieves accuracies of < |%. These cases 
arc sufficiently simple that they could be employed in hand calculations to 
yield reasonable estimates. We therefore tabulate the coefficients in Tables 8 
and 9, respectively. Increasing the number of parameters obviously improves 
the accuracy and we would judge that the M = 3 = 4 case is sufficient 
for most purposes. Here the MM deviations are < 0.02% and the RMS 
deviations arc ~ 0.003%. The 15 + 1 parameters arc given in Table 10. 

As in the case of fermions, it might be preferable, or even necessary, to 
have an approximation which merges exactly into the asymptotic limiting 
cases. Thus, it is natural to employ the values of Table 6 as constraints 
and fit the remaining parameters. For M — 3 and N — A this leaves 9+1 
free parameters. The deviations, denoted by a dagger in Table 7, are, of 
course, larger than in the unconstrained case, but only by a factor of < 2 
(the degradation was slightly worse for fermions). The coefficients for this 
case, tabulated in Table 11, are similar to those of Table 10 since the errors 
do not differ greatly. 

5 CONCLUDING REMARKS 

We have generated polynomial approximations to the thermodynamic inte- 
grals which embody and merge exactly into the analytic limiting cases and 
accurately interpolate over the entire range of temperature and chemical 
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potential. The boson case had not been treated previously and we have im- 
proved upon the parameterization of EFF and the modification of Pols et 
al. (1995) for fermions. With 10 (11) free parameters, and 20 (16) terms, 
we achieve accuracies of better than 0.008% (0.03%) for the boson (fermion) 
pressure. Some simpler, but less accurate, approximations were also dis- 
cussed. 

In astrophysics, it is often the case that particles and antiparticles are in 
thermal equilibrium. Then, the net number density is A/" = A/+ — A/l where 
Af+{Af-) is the number density of particles (antiparticles), and the total pres- 
sure is P = P+ + P_. In thermal equilibrium, we must have /i+ = — The 
problem usually posed is such that one must determine the thermodynamic 
quantities given Af and T. An obvious approach would be to employ the EFF 
scheme and to establish the limits of the pressure in analogy to equation (|^) 
in the fermion case. Unfortunately, the non-degenerate, non-relativistic case 
has no single limit in terms of power laws because of pair formation. This is 
true in both the fermion and boson cases. Thus, the EFF scheme cannot be 
directly applied to pairs in equilibrium. 

It is necessary, therefore, to treat pairs in terms of the particles and an- 
tiparticles separately. Thus, we can use the results of §3 for the number den- 
sities A/'+(/+,T) and A/1(/_,T) of fermions and antifermions, respectively. 
Equation ([T5|) gives ip± in terms of f±. The pair of simultaneous equations 

Ar = Ar+(/+,T)-Ar_(/_,T) ; ^+(/+) = -^_(/_)-2mcVT,(29) 

is then solved for /+ and as opposed to the inversion of a single equation 
M = J\f{f, T) for /. For bosons, one simply replaces f± by h± in equation 
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Table 1. Fermion coefficients, pmn, derived from the leading 
terms in the limiting expressions for p and s. 



Pmn 


n = 




n = N 


m — 


e2 nr 

a V 32 




2a 


m = M-l ^a~3 


[Stt^o + 8(4M - 5a - 1)] 




-4) + 2(M-l)] 


m = M 






4 

3a2 



Table 2. Maximum modulus (MM) and root mean square (RMS) 
deviations for the fermion thermodynamic quantities for various 
dimensionalities M and A^. 



a fitted a = 1 



p u p p 



M 


N 


MM 


RMS 


MM 


RMS 


MM 


RMS 


MM 


RMS 


2 


1 


2.5(-2) 


9.5(-3) 


3.1(-2) 


1.0(-2) 


2.5(-2) 


9.1(-3) 


2.5(-2) 


1.0(-2) 


2 


2 


8.1(-4) 


3.4(-4) 


8.0(-4) 


3.3(-4) 


9.1(-4) 


3.6(-4) 


2.1(-3) 


7.6(-4) 


3 


2 


5.9(-4) 


2.1(-4) 


5.2(-4) 


1.9(-4) 


5.9(-4) 


2.0(-4) 


7.4(-4) 


3.1(-4) 


2 


3 


7.0(-4) 


2.7(-4) 


7.0(-4) 


2.6(-4) 


9.0(-4) 


3.0(-4) 


1.9(-3) 


7.1(-4) 


3 


3 


7.1(-5) 


2.6(-5) 


1.3(-4) 


4.7(-5) 


1.3(-4) 


4.0(-5) 


5.6(-4) 


1.8(-4) 


4 


3 


7.2(-5) 


2.3(-5) 


1.3(-4) 


4.6(-5) 


7.2(-5) 


2.4(-5) 


1.9(-4) 


7.4(-5) 


3 


4 


4.0(-5) 


1.2(-5) 


3.8(-5) 


1.2(-5) 


1.2(-4) 


3.3(-5) 


5.4(-4) 


1.8(-4) 


t3 


3 


3.0(-4) 


8.1(-5) 


4.5(-4) 


1.0(-4) 


3.0(-4) 


8.7(-5) 


2.6(-3) 


8.6(-4) 



The notation 5.9(— 4), for example, indicates 5.9 x 10 ^. 
t Employing the constraints of Table 1. 
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Table 3. Fermion coefficients Pmn for M = = 2; a = 0.442. 



Pmn 




n = 


n = 1 


n = 2 


m — 





5.23810 


12.4991 


8.36058 


m — 


1 


11.1751 


25.3687 


15.6453 


m — 


2 


5.91800 


12.4945 


6.82530 



Table 4. Fermion coefficients ioi M — N — 3; a — 0.420. 



Pmn 




n = 


n= 1 


n = 2 


n — 3 


m — 





5.51219 


18.6215 


22.0078 


8.71963 


m — 


1 


17.4343 


57.7188 


66.1098 


25.6323 


m = 


2 


18.1239 


58.7781 


65.1429 


24.4246 


m = 


3 


6.30952 


19.8967 


21.1375 


7.55866 



Table 5. Fermion coefficients Pmn ior M — N — 3 employing the 
constraints of Table 1; a = 0.433. 



Pmn 




n = 


n=l 


n = 2 


n = 3 


m = 





5.34689 


18.0517 


21.3422 


8.53240 


m — 


1 


16.8441 


55.7051 


63.6901 


24.6213 


ni = 


2 


17. 1708 


56.3902 


62.1319 


23.2602 


m — 


3 


6.07364 


18.9992 


20.0285 


7.11153 
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Table 6. Boson coefficients, Pmn, derived from the leading terms 
in the hmiting expressions for P and p. 



Pmn 


n = 




III = 




2C(1) 


m — 1 


^/f(^ + l)C (1) 


2(M+1)C(4) 


m = 2 


-a-iC(i)+M(M + l)C(|); 


-a-iC(3) + M(M + l)C(4) 




p 

2 


2e^/a 



Table 7. Maximum modulus (MM) and root mean square (RMS) 
deviations for the boson thermodynamic quantities for various 
dimensionahties M and N. 



pup 



M 


N 


MM 


RMS 


MM 


RMS 


MM 


RMS 


2 


1 


2.1(-2) 


1.3(-2) 


2.3(-2) 


1.3(-2) 


4.4(-2) 


1.5(-2) 


2 


2 


2.7(-3) 


l.l(-3) 


2.9(-3) 


1.2(-3) 


5.7(-3) 


1.7(-3) 


3 


2 


3.9(-4) 


1.2(-4) 


4.6(-4) 


1.9(-4) 


3.0(-3) 


7.5(-4) 


2 


3 


2.6(-3) 


7.9(-4) 


2.7(-3) 


1.0(-3) 


4.9(-3) 


1.4(-3) 


3 


3 


l.l(-4) 


5.8(-5) 


2.0(-4) 


8.4(-5) 


5.7(-4) 


1.6(-4) 


4 


3 


l.l(-4) 


3.7(-5) 


2.4(-4) 


5.8(-5) 


5.9(-4) 


1.5(-4) 


3 


4 


4.9(-5) 


1.7(-5) 


1.3(-4) 


3.5(-5) 


1.7(-4) 


4.7(-5) 


t3 


4 


7.8(-5) 


2.4(-5) 


1.7(-4) 


4.3(-5) 


2.4(-4) 


9.0(-5) 



t Employing the constraints of Table 6. 



22 

Table 8. Boson coefficients for M = 2 and N = 1; a = 0.978. 



Prnn 




n = 


n = 1 


m = 





1.63146 


2.11571 


m — 


1 


Pin = 


3pon 


m — 


2 


3.31275 


5.15372 



Table 9. Boson coefficients Pmn for M = A?" = 2; a = 0.914. 



Pmn 


n = 


n= 1 


n = 2 


m — 





1.68131 


3.47558 


2.16582 


m = 


1 




Pin = 3pon 




m = 


2 


3.25053 


7.82859 


5.19126 



Table 10. Boson coefficients p^n for M = 3, iV = 4; a = 1.029. 



Pmn 


n = 


n= 1 


n = 2 


n = 3 


n = 4 


m — 





1.68134 


6.85070 


10.8537 


7.81843 


2.16461 


m — 


1 






Pin = 4pon 






m = 


2 


8.19651 


35.6058 


57.7131 


12.3593 


11.8199 


m — 


3 


3.45614 


15.1152 


25.5254 


19.2745 


5.51757 
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Table 11. Boson coefficients pmn for M = 3, N = A including 
the constraints of Table 6; a = 1.040. 



Pmn 




n = 


n — 1 


n = 2 


n — 3 


n = 4 


m — 





1.68130 


6.85060 


10.8539 


7.81762 


2.16465 


m — 


1 






Pin = 4pon 






m — 


2 


8.51373 


35.6576 


57.7975 


42.4049 


11.8321 


m — 


3 


3.47433 


15.1995 


25.6536 


19.3811 


5.54423 
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